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Abstract 


A method is described for evaluating the minimum distance parameters of 
trellis-phase codes, including CPFSK, partial-response FM and, more 
importantly, coded CPM (continuous-phase modulation) schemes. The algorithm 
provides dramatically faster execution times and lesser memory requirements 
than previous algorithms. Results of sample calculations and timing 
comparisons are included. 


AN IMPROVED ALGORITHM FOR EVALUATING TRELLIS PHASE CODES 


I . INTRODUCTION 

Trellis phase codes a?:e digital phase-modulated (constant-envelope) 
signals whose phase trajectories undergo smooth phase transitions which can 
have memory induced by partial-response modulation and/or convolutional 
pre-coding, [1], [2]. This combination of smoothness and memory provides 
constant-envelope designs which have intrinsically good power spectra and 
enhanced energy efficiency. The terra "trellis" arises from the finite-state 
trellis description which is possible for rational modulation indices, as 
would be the case in any practical implementation. 

Study of such signals for use on the white Gaussian noise channel with 
coherent detection is based on finding the minimum distance, in the L 2 ’“norm, 
between any pair of signals Sj(t) and S 2 (t). This is a considerably more 
difficult problem than finding the minimum Hamming distance of a convolutional 
code since the group property is Lacking and each "transmitted" sequence may 
have a unique set of distances to all "received" sequences. Thus to 
exhaustively determine minimum distance for signals of length N symbols with 
M-ary modulation requires examination of M^^ pairs of sequences of length N, 
or calculations of distance increments, which can be viewed as a basic 

unit of computation. 

If one considers a phase tree in which the phase trajectories of all 
signals are plotted as a function of time, it is possible to imagine two 
signals having the same trajectories up to a certain point, their paths 
diverging then meeting again at some later time. It is desired to determine 
the minimum distance between any pair of signals which have split and 


reme rged. 
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Previous approaches to finding the minimuin distance have in fact explored 
all pairs of paths in the phase tree which split at the initial tree level . 
Pruning rules have been added so that if the cumulative distance exceeds the 
present minimum distance for any merged pair, then descendants of the current 
pair need not be considered. Pairs which are still contenders for 
establishing the free distance, defined as the absolute minimum distance for 
all pairs of signals, are kept alive by putting them on a memory stack l3] or 
by use of a forward/backward search in the tree, [4J . The former is fast 
relative to the latter, but can be quite consumptive of memory. Neither, 
however, exploits the important finite-state nature of the signal (for 
rational modulation indices), which is the intent of this paper. 

Specifically, we show an algorithm for finding the distance to depth N 
whose complexity is upper-bounded by ( 5 ( 2 ^^ • • q • N) units of computation, 

which is linear in N although the multiplier constant may be large. In the 
above 2'^ is the number of "data states" in the encoder/modulator, M is the 
number of signalling options in each interval, and q is the denominator of the 
modulation index, h * p/q* Furthermore, memory resources necessary are very 
manageable, <5(22'' • q) in contrast to the potentially large stack size 
required for previous tree-searching algorithms, [3], 

We shall describe and illustrate the algorithm in the context of 
convolutionally-coded CPFSK, [2], [4l , although we note that the algorithm may 
be applied Co any finite-state phase-coded modulation scheme. These range 
from MSK on the simple side to M-ary partial-response FM [1] and mulci-h codes 
[5l on Che more complicated side. A necessary restriction is that h be 
rational. Obviously any irrational h can be closely approximated by a 
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rational h, with sufficiently large q, making the trellis size large. We 
don't feel this is any practical limitation, however, for actual 
implementations of modems will be restricted to rational h with small p and q 
for the same complexity reasons. 


ORIGINAL PAGE IS 
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11. BACKGROUND 

Let the modulated signal over the interval nT < t < (n+l)T be described as 

S(t, ^ “ ( 2E ^ 1 / 2 jjQg (u)(,t + (}>(t, b)) (1) 

T 

where ij)(t, b) is the phase modulation in response to the data sequence b “ (bj, 
b 2 > ••• bjj). E is the energy per symbol, and T is the symbol interval. We may 
write the phase in general as 

({>(t, b) = ^n-1 * ^n-v+l» ®n^ ™ ^^^n> ®n* ^n^ ^2) 

where is the "data state" comprised of a finite number of previous data 
symbols and "phase state", which amounts to the cumulative phase 

induced by all data symbols whose influence has ceased (see e.g. [1], [4]), 

In the unmodulated CPFSK case for example, 

n-1 (3) 

0^ = 2 b* hTT 

j=l 

The "state" of the modulator is defined as = (S^, 9^), and the structure 
of the coding/modulation scheme induces a known state-transition equation. For 
h = p/q, a finite-state representation having * q states follows, with v 
playing the role of "constraint length", or memory length of the combined 
coder/modulator. 

To decode such a signal in maximum likelihood fashion, the Viterbi 
algorithm may be applied to a trellis having 2'^ • q states. For the Gaussian 
channel, path metrics are correlations between the received signal and the 
hypothesized transmitted signal. Our interest however is in determining the 
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asymptotic performance of such a receiver, which for the assumed channel 
requires determining the minimum distance between any two signals which are 
merged at some time, split, then reraerge later. Or, for finite-memory 
receivers, we seek the minimum distance between pairs which split and are of 
length N, denoted djj. We also define the free distance 

«lfree “ 

N-^09 


as usual. 

The appropriate distance measure for this problem is 

NT 

d2 (N) = / [Si(t) - S2(t)]2dt (5) 

0 

which reduces to, for >> Zir/T, 

NT 

d2 (n)= ( 2E^ J [1 _ cos A(j>(t)]dt (6) 

12 TO 

where A<j)(t) is the time-varying phase separation between the two signals of 
length NT seconds. Observe that the value for df^gg depends on h, the 
modulation index; <})(t,b), the phase modulation charactersitic ; and the 
"memory" of the encoding process which translates into longer remerger times. 
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HI. ALGORITHK DESCRIPTION 

In order to examine all pairs of sequences we form a pair-state trellis. 

Each node in the trellis is defined by a triple consisting of Sj and S 2 > 

n n 

the two data states, and * ^1 ~ ^ 7 * phase difference at the start 

n n 

of the current interval. It is only necessary to specify the modulo 2ir phase 
difference which in fact allows us to use a finite-state model and a 
trellis. It may be readily shown that there are q such phase differences, 
where again h =■ p/q. Thus if there are 2'^ data states we now have q . 2^ 
pair-states in the trellis. The transitions from states at level n to 

states at level n+1 are defined by the encoder/ modulator structure and have 
incremental distances associated with them. 

We note that the cumulative squared distance may be recursively computed 
as 

(n+l)T 

d^l 2 (n-fl) = (n) + (^) / [1 - cos Ai^(t)]dt (7) 

T nT 

and is completely specified by transitions in the pair-state trellis. The 
task now is to find the smallest distance between any pair-state corresponding 
to a merger at level n = 0 and another merged pair-state at some level n, 
while transitioning between pair-states in the pair-state trelLis. This is 
completely analogous to the shortest-route problem and may be solved by 
application of the dynamic programming solution in the form of the Viterbi 
algorithm, [6] . The "principle of optimality" here is that if a sequence pair 
(bj, b 2 ) is to produce the minimum distance event, it wilt do so via 
extensions of minimum distance pairs to some intermediate pair-state, and this 
holds for all n. Thus it is sufficient to preserve the information associated 
with the minimum distance pair for each pair state at each level n, and 
proceed forward recursively using the known pair-state transitions. 
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To locate Che mlnlioum distance for each pair state j we fetch previous 
distances for the pair states which transition to j, add Che incremental 
distances associated with chose transitions, make roughly comparisons to 
find the minimum, and save this distance. If it is desired to locate the 
minimum-producing sequence pair, path maps can be stored and updated as well. 
Since there are 2 ^''^ . q pair-states and branches/s Cate we say the 
conq>lexity is • q • • N) for an N-level trellis. 

Actually the same pruning rules used in earlier implementations may be 
applied Co speed Che calculation. Specifically, some initial (optimistic) 
estimate of d^fj-ee is entered, and if the minimum distance of a pair leading 
Co a pair-state merge in the trellis (of which there are 2 ^'^ at every level) 
falls below the initial estimate, the d^f^.g,g value is decreased to the newly 
found value, d^fj-gg is decreased further if other subsequent merges give 
lesser distances. Simultaneously, Khenever the cumulative distance to a 
non-merged pair-state exceeds the current this pair-state need not be 

extended in the next trellis level, for descendants of this pair~state cannot 
ultimately produce the event. 

The distance profile as a function of n, d2j^i„(n), is also easily 
obtained by locating the minimum of distances over all survivor pairs to 
states at level n. 

To initialize the algorithm, we assign zero distance Co all merged 
pair-states in the trellis, and large initial distance Co all others. In 
effect this allows the search for pair sequences beginning with any merged 
state, and thus Che examination of all possible starting states is 
incorporated automatically. In contrast, the tree-searching procedures have 
simply rooted the tree search in each possible initial state and found a 
global distance minimum over this set of states. 
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A flow-chart of the algorithm is shown in Figure 2. To give better feel 
for the level of complexity, we list in Table 1 the trellis parameters for a 
variety of pertinent cases. Application of the algorithm to seemingly 
disparate cases is rather easy: all that must be done is identify a proper 

state-description and its associated state transitions, and to define a 
case-specific subroutine which computes a table of incremental distances for 
various pair-state transitions. 

The comparison between the trellis-based algorithm described here and the 
tree-searching algorithms used previously is not as one-sided (exponential 
versus linear in N) as the above complexity relations indicate. Pruning rules 
in both algorithms speed the execution substantially, and in fact, [1] reports 
search time much faster than exponential in N. Still the trellis algorithm 
will be faster since it exploits pair-state merges while the tree-based 
algorithms do not. In addition, our memory requirements are generally much 
smaller, being known in advance as well, in contrast to the less predictable 
stack-size requirement of tree-searching methods. Specific comparisons are 
offered in the next section. 
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IV. EXAMPLES AND DISCUSSION 

Out present interest is in evaluating the use of convolutional coding 
combined with CPFSK and other CPM methods as a combined power/bandwidth 
efficient scheme. The large number of possible code and modulation paramet<irs 
makes an efficient distance computation essential. 

First we consider the example of the rate 1/2, constraint-length 4 (v » 

3) encoder of Figure 3. The two output code symbols are treated as a single 
quaternary symbol in uhe set {-3, -1, 1, 3) with natural-binary mapping 
assumed, i.'e. 00 - -3, 11 *• +3. We have earlier found this scheme to be 
optimal in its class for h < 1/4. The free distance as a function of h for 
this code is plotted in Figure 3, where we see monotone-increasing 
to h slightly beyond 1/4. Thereafter other merger events dominate the 
distance, and this code is no longer optimal, [4] . For certain values of h, 
1/4 and 1/2 in this case, the distance has isolated small values owing to 
unusually short merger events. These are known as weak modulation indices, 
and other codes perform better at these h. 

Also shown in Figure 4 and 5 are the distance results for r “ 2/3, v * 2 
coded octal CPFSK (the optimal small h code in this class [7]) and for 4-ary 
uncoded CPFSK, for which the distance result was already known. Both schemes 
fit into the trellis phase code framework and are also readily evaluated by 
this procedure. These distance plots were obtained by computing distance for 
all rational h with h < l/2 and q < 50. Due to this discretization, the 
distance near the weak modulation indices is not represented exactly. 

Also readily available from the algorithm is the "distance profile," 
djnin^(n) versus n for various h, or vice-versa. This reveals how much 
receiver memory is necessary so that all unraerged pairs have distance 
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exceeding For example, for the r * 1/2 code, with h < 1/4, 11 symbol 

delay is sufficient, while for r * 2/3, h < 1/8, 12 symbol delay is adequate. 

Table 11 lists execution times for the algorithm described here and a 
tree'-searching procedure similar to that of [3] . Both were programmed in 
FORTRAN by reasonably competent programmers, using similar data types and 
structures. We see the new method is faster as predicted, although the tree 
method's pruning rules allow it to be reasonably close. The speed advantage 
diminishes as q increases due to the increasing trellis size. In this case, 
many of the pair-states are "dead", but merely checking this fact reduces the 
efficiency somewhat. As a general conclusion we can say the trellis algorithm 
has its biggest gain when the tree is in fact mergeable into a small trellis, 
which occurs for smaller memory lengths and q small. 
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Signal Type # pair-states branches/state 


MSK, h = 1/2 

2 

(2) 

4 

4-ary CPFSK, h = 1/4 
(uncoded) 

4 

(4) 

16 

4-ary 2 RC, h = 1/4 

64 

(40) 

16 

8-ary CPFSK, h = 1/8 

8 

(8) 

64 

r » 1/2, V * 2, coded 
4-ary CPFSK, h » 1/4 

64 

(40) 

4 

r = 2/3, V * 2, coded 
8-ary CPFSK, h = 1/8 

128 

(80) 

16 

binary {4/8, 5/8} raulti-h 

8 

(8) 

4 


code 


Table I. 

Pair-State Trellis Parameters for Several Trellis Phase Codes 


Case 

r 

h 

V 

Trellis Algorithm* 

Tree Algorithm* 

d^f ,.„p/2Ef^ 

I 

1/2 

1/4 

3 

13 sec. 

47 sec. 

5.24 

II 

II 

1/5 

3 

16 

32 

4.30 

III 

II 

3/16 

3 

27 

39 

3.50 

IV 

II 

1/16 

3 

21 

36 

0.48 

V 

2/3 

1/8 

2 

21 

— 

2.56 

* Run 

in FORTRAN on 

POP 

11/03 microcomputer. 




Table II. 


Comparison of Execution Times for Tree and Trellis Algorithms 
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Finite-State Machine 


Figure 1. Conceptual Description of Distance-Finding Problem: Vector- 
Input Sequence Produces Scalar-Output Distance According to 
Finite-State (Trellis) Description. 
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Figure 4. Distance versus h for r = 2/3, v * 2, coded octal CPFSK 





